Towards singularity and ghost free theories of gravity 
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We present the most general covariant ghost-free gravitational action in a Minkowski vacuum. 
Apart from the much studied f(R) models, this includes a large class of non-local actions with 
improved UV behavior, which nevertheless recover Einstein's general relativity in the IR. 



The theory of General Relativity (GR) has an ultravi- 
olet (UV) problem which is typically manifested in cos- 
mological or black-hole type singularities. Any resolution 
to this problem requires a theory which is well behaved 
in the UV and reduces suitably to Einstein's gravity in 
the infrared (IR) 1 . In this letter, our aim is to investigate 
whether the typical divergences at short distances can be 
ameliorated in higher derivative covariant generalizations 
of GR. 

Higher derivative theories of gravity arc generally bet- 
ter behaved in the UV and offer an improved chance 
to construct a singularity free theory Furthermore, 
Ref. [H demonstrated that fourth order theories of grav- 
ity are renormalizable, but inevitably suffer from unphys- 
ical ghost states. Therefore, before we address the short- 
distance behavior of GR, we first ennumerate the subset 
of all possible modifications to Einstein's gravity which 
are guaranteed to be ghost-free. To the best of our knowl- 
edge, a systematic method for this is not presently avail- 
able. 

Generic quadratic action of gravity: Let us start 
with the most general covariant action of gravity. We im- 
mediately realize that to understand both the asymptotic 
behavior in the UV and the issue of ghosts, we require 
only the graviton propagator. In other words, we look at 
metric fluctuations around the Minkowski background 



(1) 



and consider terms in the action that are quadratic in 
h^. Since in the Minkowski background R^ v \ a vanishes, 
every appearance of the Riemann tensor contributes an 
0{h) term in the action. Hence, we consider only terms 
that are products of at most two curvature terms, and 
higher ones simply do not play any role in this analysis. 



1 In the light of current cosmic acceleration observations, there 
have been efforts to modify gravity at large distances, see pj for 
but we do not discuss these models here. 



The most general relevant action is of the form 
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x v ~9 
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where O is a differential operator containing covariant 
derivatives and r/^. We note that if there is a differ- 
ential operator acting on the left Riemann tensor, one 
can always recast that into the above form by integrat- 
ing by parts. The most general action is captured by 
14 arbitrary functions, the F^s, which reduce to the 6 
we display in eq. (|27|) upon repeated application of the 
Bianchi identities. 

Our next task is to obtain the quadratic (in h^ v ) free 
part of this action. Since the curvature vanishes on the 
Minkowski background, the two h dependent terms must 
come from the two curvature terms present. This means 
the covariant derivatives take on their Minkowski values. 
As is obvious, many of the terms simplify and combine 
to eventually produce the following action 
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h <T b{U)d a d v hf iV (3) 



hc{U)d ll d v h> lv + -hd(D)ah 
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The above can be thought of as a higher derivative gener- 
alization of the action considered by van Nieuwenhuizen 
in Ref. Here, we have allowed a,b,c,d and / to be 
nonlinear functions of the derivative operators that re- 
duce in the appropriate limit to the constants a, 6, c and 
d of Ref. Q. The function /(□) appears only in higher 
derivative theories. In the appendix (j29H33[) we have cal- 
culated the contribution from the Einstcin-Hilbert term 
and the higher derivative modifications to the action in 
eq.Q. From the explicit expressions we observe the fol- 
lowing relationships: 



a + b = 
c + d = 
b+c+f = 



(4) 
(G) 



a review. 



so that we are left with only two independent arbitrary 
functions. 
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The field equations can be derived straightforwardly to 
yield 

a(npv + b(a)d a (d^ l + d^) 

+ c{a)( v ^d p d a hr° + dpdvh) + v^d(a)ah 
+ fi^n-^dxd^h^ = -KT^ . (7) 

While the matter sector obeys stress energy conservation, 
the geometric part is also conserved as a consequence of 
the generalized Bianchi identities: 



ktV, 



r£ = = (a + b)DK, 



M + (c + d)Ud„h 
f)K&v (8) 



It is now clear why eqs.(@][l)]) had to be satisfied. 
Propagator and physical poles: We are now well- 
equipped to calculate the propagator. The above field 
equations can be written in the form 



(9) 



where 11^ is the inverse propagator. One ob- 
tains the propagator using the spin projection operators 
{P 2 ,P S °,P°,P^}, see Ref. 0. They correspond to the 
spin- 2, the two scalars, and the vector projections, re- 
spectively. These form a complete basis. Considering 
each sector separately and taking into account the con- 
straints in cq.(|4][6]), we eventually arrive at a rather sim- 
ple result 
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We note that the vector multiplet and the w-scalar have 
disappeared, and the remaining s-scalar has decoupled 
from the tensorial structure. Further, since we want to 
recover GR in the IR, we must have 



o(0) = c(0) = -6(0) = -d(0) = 1 



(11) 



corresponding to the GR values. This also means that as 
k 2 — > we have only the physical graviton propagator: 



lim IP" 



Act 



(p 2 A 2 ) - 0P s °/2fc 2 ) 



(12) 



A few remarks are now in order: First, let us point out 
that although the P s residue at k 2 = is negative, it is 
a benign ghost. In fact, P s ° has precisely the coefficient 
to cancel the unphysical longitudinal degrees of freedom 
in the spin two part Thus, we conclude that pro- 
vided cq.lfTTj) is satisfied, the k 2 = pole just describes 
the physical graviton state. Secondly, cq. pTj) essentially 
means that a and c are non-singular analytic functions 
at k 2 = 0, and therefore cannot contain non-local inverse 
derivative operators (such as a(D) ~ 1/D). 

Let us next scrutinize some of the well known special 
cases: 



/(P) gravity: they are a subclass of scalar-tensor theo- 
ries and are studied in great detail both in the context of 
early universe cosmology and dark energy phenomenol- 
ogy. Here, only the Pi appears as a higher derivative 
contribution (see appendix). According to our preced- 
ing arguments, we obtain the physical states from the 
P 2 term. Since a = 1, it is easy to see that only the s- 
multiplct propagator is modified. It now has two poles: 

II l/2fc 2 (fc 2 - m 2 ) + .... The k 2 = pole has, as 

usual, the wrong sign of the residue, while the second pole 
has the correct sign. This represents an additional scalar 
degree of freedom confirming the well known fact 0, @] . 
Fourth order modification in R pu R^ v : They have 
also been considered in the literature. This corresponds 
to having an P2 term (see appendix) , which modifies the 
spin-2 propagator: II ~ P 2 /fc 2 (fc 2 — to 2 ) + . . . . The sec- 
ond pole necessarily has the wrong residue sign and cor- 
responds to the well known Weyl ghost, Refs. 0, @]. In 
fact, this situation is quite typical: /(P) type models 
can be ghost-free, but they do not improve UV behavior, 
while modifications involving P^Ao-'s can improve the 
UV behavior 0] but typically contain the Weyl ghost! 

To reconcile the two problems we now propose first to 
look at a special class of non-local models with / = or 
cquivalently a = c. The propagator then simplifies to: 



IT"' i - - I P 2 - - P f 

Xa ~ k 2 a(-k 2 ) { 2 s 
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It is obvious that we are left with only a single arbitrary 
function a(D), since now a = c = —b = —d. Most impor- 
tantly, we now realize that as long as a(D) has no zeroes, 
these theories contain no new states as compared to GR, 
and only modify the graviton propagator. In particular, 
by choosing a(D) to be a suitable entire function we can 
indeed improve the UV behavior of gravitons without in- 
troducing ghosts. This will be discussed below. 
Singularity free gravity: We now analyze the scalar 
potentials in these non-local theories, focussing partic- 
ularly on the short distance behavior. As is usual, we 
solve the linearized modified Einstein's equations (0 for 
a point source: 



(14) 



Next, we compute the two potentials, ^(r), \&(r), corre- 
sponding to the metric 



ds 2 = -(1 + 2$)d< 2 + (1 - 2^)dx 2 



(15) 



Due to the Bianchi identities 0, H[ , we only need to solve 
the trace and the 00 component of eq. ([7]) . Since the New- 
tonian potentials are static, the trace and 00 equation 
simplifies considerably to yield 

(a - 3c)Dft + (4c - 2a + f)d lt d v h! lv = up 

aUh M + cnh-cd tt d v hf lv = -up, (16) 
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which for the metric eq. (|15p simplify to 

2(a-3c)[V 2 <J>-4V 2 *] = up 
2(c-a)V 2 $-4cV 2 * = -Kp. 



(17) 



We are seeking functions c(D) and a(D), such that there 
are no ghosts and no 1/r divergence at short distances. 
For / = 0, the Newtonian potentials are solved easily: 

4a(V 2 )V 2 $ = 4a(V 2 )V 2 * = up = nmS 3 (f). (18) 

Now, we know that in order to avoid the problem of 
ghosts, a(D) must be an entire function. Let us first il- 
lustrate the resolution of singularities by considering the 
following functional dependence 0] : 



o(D) 



-a/iu 2 



(19) 



Such exponential kinetic operators appear frequently in 
string theory In fact, quantum loops in such stringy 
non-local scalar theories remain finite giving rise to in- 
teresting physics, such as linear Regge trajectories fioj 
and thermal duality [lljj. We note that there are a wide 
range of allowed possible energy scales for M, including 
roughly the range between A and M p i. 

Taking the Fourier components of cq. (1181) , in a straight 
forward manner one obtains 



in 



Airm f dp sin pi 



(b( r ) r~j "" / rf 3 w 
{1 M| J P p 2 a(-p 2 ) rM 2 J p a(-p 2 )' 

(20) 

We note that the 1/r divergent piece comes from the 
usual GR action, but now it is ameliorated. For eq. (JTUJ) 
we have 



$(r) 



M 2 r 



P 



sin (p r) 



2M* r 



erf 



rM 
~2 
(21) 

and the same for \&(r). We observe that as r — > oo, 
erf (r) — > 1, and we recover the GR limit. On the other 
hand, as r — > 0, erf(r) — > r, making the Newtonian 
potential converge to a constant ~ mM/M 2 . Thus, 
although the matter source has a delta function singu- 
larity, the Newtonian potentials remain finite! Further, 
provided mM <C M p , our linear approximation can be 
trusted all the way to i — > 0. 

Let us next verify the absence of singularities in the 
spin-2 sector. This will allow us, for example, to derive 
a singularity free quadrupole potential. We enforce the 
Lorentz gauge as usual so that the generalized field equa- 
tions (J7J read 



lUh 



f c 

-d^dyh - -rj^Dh = -KTpv . (22) 



Again for / = we have a simple wave equation for 
the graviton a(D)D/i /il , = —kt^ v . We invert Einstein's 
equations for to obtain the Greens function, G^, for 



a point-like energy-momentum source. In other words, 
we solve for 

a(D)DG^(x ~y) = -kt„ v S\x - y), (23) 
Under the assumption of slowly varying sources, one has 

'rM' 



Gn V {r) — -7rerf 



(24) 



for a(D) given in ea. p^|) . We observe that in the limit 
r — > 0, the Greens function remains singularity free. The 
improved scaling takes effect roughly only for r < 1/M. 
Cosmological Singularities: The very general frame- 
work of this paper allows us to consistently address the 
singularities in early universe cosmology. As an example, 
we note that a solution to eq.([7]) with 

h ~ diag(0, A sin At, A sin Xt, A sin Xt) with A < 1 (25) 

describes a Minkowski space-time with small oscillations 
fl2j . This configuration is singularity free. Evaluat- 
ing the field equations for cq. (|2"5)) gives the constraint 
o(— A 2 ) — 3c(— A 2 ) = 0. Thus, our simple / = case 
is not sufficient and we require an additional scalar de- 
gree of freedom in the s-multiplet. Note that this also 
explains why a solution such as cq. (|25[) is absent in GR. 
We generalize to / ^ 0, but take special care to keep in- 
tact our results in eq. (fTTj) and eq. (|18p . The most general 
ghost-free parameterization for a ^ c is 



c(D) 



_ «(□) 



1 + 2 1 



(26) 



where c(D),a(D) are entire functions. Note that m 2 — > 
oo and c = 1 reproduces the / = limit. We now find 
that eq. ([25[) is a solution to the vacuum field equations 
with A = m. How the universe can grow in such models 
and also how the matter sector can influence the dynam- 
ics can possibly be addressed only with knowledge of the 
full curvature terms. We hope to investigate this in future 
work, but see Ref. [13j and 141 for similar considerations. 



Generality: How general are the above arguments lead- 
ing to a lack of singularities? According to the Weier- 
strass theorem any entire function is written as a(0) = 
g-7(n) j wnere 'y(D) is an analytic function. For a polyno- 
mial 7(D) it is now easy to see that if 7 > as □ — > 00, 
the propagator is even more convergent than the expo- 
nential case leading to non-singular UV behavior. 
Conclusion: We have shown that by allowing higher 
derivative non-local operators, we may be able to render 
gravity singularity free without introducing ghosts or any 
other pathologies around the Minkowski background. It 
should be reasonably straight-forward to extend the anal- 
ysis to DcSittcr backgrounds by including appropriate 
cosmological constants. In fact, requiring that the the- 
ory remains free from ghosts around different classical 
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vacua may be a way to constrain the higher curvature 
terms that didn't seem to play any role in our analysis. 
Other ways of constraining/determining the higher cur- 
vature terms would be to look for additional symmetries 
or to try to extend Stclle's renormalizability arguments 
to these non-local theories. Efforts in this direction have 
been made (l5| . Finally, it is known that one can ob- 
tain GR starting from the free quadratic theory for h^ v 



by consistently coupling to its own stress energy tensor. 
Similarly, can one obtain unique consistent covariant ex- 
tensions of the higher derivative quadratic actions that 
we have considered? We leave these questions for future 
investigations. 
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itational action. 



Appendix 

The quadratic action in curvature reads 

S q = J d^x^lRF^R + R^F 2 (a)R^ + R^ Xa F 3 (n)R^ Xtr + i?F 4 (D)V Al V I ,V A V (T .R^ 

+ i?^ lCTl F 5 (D)Vp 1 V CTl V, 1 V,V p V CT i?'" yACT + iF 1 ^ lCTl F 6 (D)V Pl V CTl V, 1 V Pl V p V,VpV ff i?^ Aff ] (27) 
where we have used the Bianchi identities: 

^crR^vXp + V ' pRpvaX + V \R/j, U pcT = (28) 

to absorb all the other covariant terms into the above six. Further, in the F4, F 5 and F 6 terms, one ends up with 
anticommutator of the covariant derivatives due to the anti-symmetric properties of the Reimann tensor, but these 
anticommutators produce a third curvature term, and therefore these terms are at least 0(h 3 ). Thus, the coefficients 
of the free theory (|3|) in terms of the F's are given by 

a(D) = l-^F 2 (D)D-2F 3 (n)D (29) 

b(n) = -i + ^F 2 (n)n + 2F 3 (n)n (30) 

c(D) = l + 2Fi(n)D + ^F 2 (D)n (31) 
d(n) = -l-2fi(np--F a (D)n (32) 
/(□) = -2Fi(Dp - F 2 (Dp - 2F 3 (Dp. (33) 
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